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Abstract. Macroscopic electroencephalographic (EEG) fields can be an explicit top- 
down neocortical mechanism that directly drives bottom-up processes that describe 
memory, attention, etc. The top-down mechanism considered are macrocolumnar 
EEG firings in neocortex, as described by a statistical mechanics of neocortical 
interactions (SMNI), developed as a magnetic vector potential A. The bottom-up 
process considered are Ca 2+ waves prominent in synaptic and extracellar processes 
that are considered to greatly influence neuronal firings. Here, the complimentary 
effects are considered, i.e., the influence of A on Ca 2+ momentum, p. The canonical 
momentum of a charged particle in an electromagnetic field, II = p + qA (SI units), is 
calculated, where the charge of Ca 2+ is q = — 2e, e is the magnitude of the charge of an 
electron. Calculations demonstrate that macroscopic EEG A can be quite influential on 
the momentum p of Ca 2+ ions, in both classical and quantum mechanics. Molecular 
scales of Ca 2+ wave dynamics are coupled with A fields developed at macroscopic 
regional scales measured by coherent neuronal firing activity measured by scalp EEG. 



PACS numbers: 87.19.L-, 05.10.Gg, 87.50.yg. 03.67.Hk 
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1. Introduction 

1.1. Multiple Scales 

There is a growing awareness of the importance of multiple scales in many physical 
and biological systems, including neuroscience (Tj, [2] . As yet, there do not seem to 
be any explicit top-down mechanisms that directly drive bottom-up processes that 
describe memory, attention, etc. Of course, there are many top-down type studies 
demonstrating that neuromodulator [3] and neuronal firing states, e.g., as defined by 
electroencephalographic (EEG) frequencies, can modify the milieu of individual synaptic 
and neuronal activity, which is still consistent with ultimate bottom-up paradigms. 
However, there is a logical difference between top-down milieu as conditioned by some 
prior external or internal conditions, and some direct top-down processes that direct 
cause bottom-up interactions specific to short-term memory (STM). 

This study crosses molecular (Ca 2+ ions), microscopic (synaptic and neuronal), 
mesoscopic (minicolumns and macrocolumns) , and macroscopic (regional scalp EEG) 
scales. Calculations support the interaction between synchronous columnar firings large 
enough to be measured by scalp EEG and molecular scales contributing to synaptic 
activity: On one hand, the influence of macroscopic scales on molecular scales is 
calculated via the evolution of Ca 2+ quantum wave functions. On the other hand, 
the influence of Ca 2+ waves is described in the context of a statistical mechanics model 
that already has been verified as calculating experimental observables, aggregating and 
scaling up from synaptic activity, to columnar neuronal firings, to regional synchronous 
activity fit to EEG while preserving an audit trail back to underlying synaptic 
interactions. 



1.2. Magnetism Influences in Living Systems 

There is a body of evidence that suggests specific top-down mechanisms for neocortical 
STM processing. 

An example of a direct physical mechanism that affects neuronal processing not 
part of "standard" sensory influences is the strong possibility of magnetic influences 
in birds at quantum levels of interaction [2}j6]. It should be noted that this is just a 
proposed mechanism (7j. 

The strengths of magnetic fields in mammalian neocortex may be at a threshold 
to directly influence synaptic interactions with astrocytes, as proposed for long-term 
memory (LTM) [8] and STM (9j[To] . Magnetic strengths associated by collective EEG 



activity at a columnar level gives rise to even stronger magnetic fields. Columnar 
excitatory and inhibitory processes largely take place in different neocortical laminae, 
providing possibilities for more specific mechanisms. 

Note that magnetic fields generated by axons, about 1CT 7 T, are generally small 
relative to the Earth's magnetic fields on the order of 3 x 1CT 5 T (T = Telsa = kg-A -1 - 
s~ 2 ). This is often cited as a reason that magnetic fields due to firing neurons are not 
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influential in brain processes. However, this paper calculates the oscillatory magnetic 
vector potential A due to many synchronous minicolumns from many macrocolumns, 
not the magnetic field of a single axon or minicolumn. The effects of this oscillatory A, 
synchronized to Ca 2+ waves which contribute to this process, are the issue. 



1.3. SMNI Context o/Ca 2+ Waves 

Since 1981, 30+ papers on a statistical mechanics of neocortical interactions (SMNI) 
applied to columnar firing states, have detailed properties of short-term memory — e.g., 
capacity (auditory 7 ± 2 and visual 4 ± 2), duration, stability , primacy versus recency 
rule, Hick's law — and other properties of neocortex by scaling up to macrocolumns 



across regions to fit EEG data [ll}|17j. SMNI not only details STM, testing SMNI at 
columnar scales, but also scaled-SMNI at relatively macroscopic scales has very well 
modeled large EEG databases, testing SMNI at regional scales. 

Experimental research supports information/memory processing by coherent 
columnar firings across many neurons [l8j[l9]. This experimental confirmation greatly 
enhances the importance the SMNI approach. There most likely are additional neural 
mechanisms that actually code information within the context of such synchronous 
firings [20] . The neuroscience community also is accepting that long-term memories 
are not stored in individual neurons, but in groups of neurons perhaps as small as 



macrocolumns 21 



The influence on the momentum of a Ca 2+ ion from macrocolumnar EEG fields as 
measured on the scalp, is thereby considered as the processing of information. SMNI 
calculates the influence of these regional synchronous firings at molecular scales that 
drive most influential Ca 2+ waves across synapses. The roles of Ca 2+ in neocortex, 
while not completely understood, are very well appreciated as being quite important 
in synaptic interactions, e.g., in modulating the production of excititory glutamic 



acid 22 , albeit not in all synpatic interactions 23 . It also is possible that Ca 



,2+ 



waves may be instrumental in tripartite synaptic interactions of astrocytes and neuronal 



synapses 24-26 , although the role of tripartite synapses in the adult brain has been 



disputed 27 . Researchers are more regularly examining glial cells to better understand 



neural processing of information 28 . The Ca 2+ waves considered here arise from a 
nonlinear cooperative regenerative process from internal stores, complementary to Ca 2+ 
released through classic endoplasmic reticulum channels and voltage-gated and ligand- 
gated Ca 2+ transients |26|. 

Columnar EEG firings calculated by SMNI lead to electromagnetic fields which can 
be described by a vector potential. This is referred to as the SMNI vector potential 
(SMNI- VP). An early discussion of SMNI- VP contained in a review of short-term 



memory as calculated by SMNI was simply suggested 17 , and a previous paper outlined 
the approach taken here, but only in a classical physics context [29]. Current research 
is directed to more detailed interactions of SMNI- VP firing states with Ca 2+ waves. 
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1.4- Outline 

Section 2 calculates the influence of A, derived using current experimental data, on Ca 2+ 
momenta p in both classical and quantum physics. Since Ca 2+ waves are influential in 
synaptic interactions, this demonstrates the influence of A in synaptic interactions. 

In Section 3 SMNI is scaled to regional EEG activity, e.g., as measured on the 
scalp, and similar to previous studies that fit synchronous columnar firings to the 
electric potential $ to EEG data, here such synchronous columnar firings are scaled 
to the vector potential A. In either case, developing $ or A, SMNI provides an audit 
trail back to columnar averaged synaptic parameters which are fit to EEG data, with 
parameters constrainted to their experimentally determined ranges. In fitting A, the 
prior established dependence of synaptic interactions on A provides justification to 
include coefficients as parameters to test and detail the dependence and senstivity of 
EEG on A. 

Section 4 discusses three current projects in the context of this paper. Some early 
results are given for one of these projects. 

The Conclusion, Section 5, summarizes the calculations. 



2. Classical and Quantum Considerations 

In these calculations, the Lagrangian formulation will be used. In descriptive terms, for 
classical physics calculations the Lagrangian L is defined by the argument of a short-time 
conditional probability distribution P over a vector of variables x and time t, 

P[x(t)\x(t - At)] = Nexp(-LAt) (1) 

where N is a normalization prefactor. This conditional probability evolves the initial 
distribution, e.g., as expressed by the path integral over all variables at all intermediate 
times. In quantum physics the Lagrangian is similarly defined in terms of the evolution 
of the wave function ip whose absolute square is a probability distribution. As 
demonstrated in many disciplines as well as in many SMNI papers, the Lagrangian 
formulation often offers intuitive, algebraic and numerical advantages to its equivalent 
partial and stochastic differential representations. E.g., this approach affords the use 
of powerful derivations based on the associated variational principle, e.g., Canonical 
Momenta and Euler-Lagrange equations. This is all rigorously derived in many 
preceding SMNI papers, and has required developing powerful numerical algorithms 
to fit these algebraic models to data, such as Adaptive Simulated Annealing (ASA) 



30 31 , to calculate numerical details of the propagating probability distributions using 



PATHINT 32,33 and PATHTREE 34. 
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2.1. Effective Momentum II 

The effective momentum, II, affecting the momentum p of a moving particle in an 



electromagnetic field, is understood from the canonical momentum 35 -38 , in SI units 



n = p + qA (2) 

where q = —2e for Ca 2+ , e is the magnitude of the charge of an electron = 1.6 x 10~ 19 C 
(Coulomb), and A is the electromagnetic vector potential. (In Gaussian units II = 
p + qA/c, where c is the speed of light.) A represents three components of a 4- 



vector [38|39 . In the standard gauge, the 3- vector components of this 4- vector potential 
related to magnetic fields are of interest. 

II can be used in quantum as well as in classical calculations. Quantum mechanical 
calculations including these effects are likely important as it is clear that in time scales 
much shorter than neuronal firings Ca 2+ wave packets spread over distances the size 
of typical synapses |40] . The gauge of A is not specified here, and this can lead to 



important effects especially at quantum scales 41 



2.2. Quantum Calculation 

The Lagrangian L, the argument of the exponential defining this probability 
distribution, includes the canonical energy II 2 /(2m). The momentum representation 
of a Gaussian wave function is developed in this context. The magentic vector potential 
field A is shown to be quite insensitive to a reasonable spatial location, so we just have 
to consider the expectation of momentum p, which essentially gives back the classical 
value. This is made more explicit as follows: 

The 3-dimensional Gaussian wave function in r-space of a Ca 2+ ion is derived as 
follows. (Further analysis is required of an initial many-bode wave function for a Ca 2+ 
wave.) The normalized wave function at time t = in momentum space for a wave 
packet centered with momentum p is 

<j>{p, 0) = (27r(Ap) 2 )- 3 / 4 e-( p - po > 2 /W Ap ) 2 ) (3) 

where squared vectors represent inner products, e.g., (Ap) 2 = Ap • Ap. There is as 
yet no experimental evidence as to how this kind of wave packet is developed by Ca 2+ 
waves in vivo. <fi develops in time as U = exp(-iHt) with Hamiltonian/Energy H, 

U(v,t) = e ^((P+9A) 2 i)/(2m/l) 

<f>(p,t) = <j>(p,0)U(p,t) (4) 

The normalized wave function in coordinate space is given by a Fourier transform in 
k-space, which can be taken in p-space using p = Kk, 

oo 

7p(r,t) = (2nh)- 3/2 [ d 3 p<j)(p,t)e ip - r/h (5) 
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This integral yields 

qAt iftpo 



3/2 







m 2(Ap) 



'itft h 2 

' 2/// 4(Ap) 2 
Pq iq 2 A 2 t 



" 4(Ap) 2 + 2mft ^ 

In coordinate space ?/> exhibits a direct dependence on A in the displacement r — > 
r — qAt/m. This is not the phase-dependence of A at work here, recognized as the 



Aharonov-Bohm effect 42 
Note that 

(Ap) 2 (Ar) 2 > (ft/2) 2 (7) 

where ((Ar) 2 ) 1//2 is the spatial 1/2-width of the packet. With the variance of if) in terms 
of 1/Ap instead of Ar, a factor of h~ 3 ^ 2 is introduced into if>(r, t) in order that the wave 
function in coordinate space be properly normalized. The dispersion of the wave packet 
in time can be seen in the factors and terms in if>(r, t), {a, (3, 7, 5}, above. 

If just the effects of A on the wave function is required, using p-space is more 
straightforward than a typical p • A calculation that does a partial integration to get 
dA/dt, giving — r • E, in terms of the coordinate r and electric field E, but r is not as 
directly observed as is p. Also note that the quantum expected value of p from J (f)*4>p 
returns just po, the same as the classical value. 

2.3. A of Wire 

For a wire/neuron carrying a current I, measured in A (not bold A) = Amperes = C/s, 

where the current is along a length z (a neuron), observed from a perpendicular distance 
r from the line of thickness tq. Neglecting far-field retardation effects, this yields 

A '"ilogf-l (9) 



47r \r 

Similar formulae for other geometries are in texts (39). The point here is the insensitive 
log dependence on distance. The estimates below assume this log factor to be of order 1. 
However, especially in this neocortical EEG context, the (oscillatory) time dependence 
of A(t) derived from I(t) is influential in the dynamics of Ca 2+ waves. 
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The magnetic field B derived from A, 

B = VxA (10) 

is still attenuated in the glial areas where Ca 2+ waves exist, and its magnitude decreases 
as inverse distance, but A derived near the minicolumns will be used there and at further 
distance since it is not so attenuated. The electrical dipole for collective minicolumnar 
EEG derived from A is 

IC If 

E = -VxB=-VxVxA (11) 

CO CO 

fj,Q, the magnetic permeability in vacuum = 47rl0~ 7 H/m (Henry/meter), where Henry 
has units of kg-m-C~ 2 , is the conversion factor from electrical to mechanical variables. 



In neocortex, jj, « 43 , 44 



The contribution to A can be viewed as including many such minicolumnar lines 
of current across 100's to 100's of macrocolumns that typically contribute to large 
synchronous bursts of EEG [45], e.g., within a region not large enough to include many 
convolutions. 

2.4- Effects of A on p 

The momentum p at issue is calculated for comparison to the vector potential. In 
neocortex, a Ca 2+ ion with mass m = 6.6 x 10~ 26 kg, has speed on the order of 
50 /im/s (46] to 100 fim/s [26]. This gives a momentum on the order of 10~ 30 kg- 
m/s. A study of molar concentrations gives an estimate of a Ca 2+ wave as comprised of 
tens of thousands of free ions representing about 1% of a released set (the bulk being 
buffered), with a range of about 100 /im, duration of about 500 ms, and concentrations 



ranging from 0.1-5 /im 26 



qA can be calculated at several scales: 



In studies of small ensembles of neurons 47 , an electric dipole moment Q is defined 
as Izr, where f is the direction unit- vector, leading to estimates of |Q| for a pyramidal 
neuron on the order of 1 pA-m = 10~ 12 A-m. Multiplying by 10 4 synchronous firings 
in a macrocolumn gives an effective dipole moment |Q| = 10~ 8 A-m. Taking z to 
be 10 2 /zm = 10~ 4 m (a couple of neocortical layers) to get I, this gives an estimate 
|gA| k2x 10- 19 x 10~ 7 x lO^/lO" 4 = 10~ 28 kg-m/s, 

Estimates at larger scales [48] give a dipole density P = 0.1 /iA/mm 2 . Multiplying 
this density by a volume of mm 2 x 10 2 /im (using the same estimate above for z), gives 
a |Q| = 10 -9 A-m. This is smaller than that above, due to this estimate including 
cancellations giving rise to scalp EEG, while the estimate above is within a macrocolumn 
(the focus of this study), leading to |gA| = 10~ 29 kg-m/s. 

The estimates for Q come from experimental data, which therefore include all 
shielding and material effects expected in other theoretical treatments that would derive 
Q. In the context of coherent activity among many macrocolumns, correlated with 



STM 19 , | A | may become orders of magnitude larger than these conservative estimates. 
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Since Ca waves play an important role in synaptic activity inherent in this coherent 
macro columnar activity, of course there is direct coherence between these waves and the 
activity of A. 

In the context of classical physics, the above calculations show how important 
the effect of qA from macroscopic EEG, on the order of 10~ 28 kg-m/s can be on the 
momentum p of a Ca 2+ ion, on the order of 10~ 30 kg-m/s. By itself, this simple numerical 
comparsion shows the important influence of A on p at classical scales. 

In the context of quantum physics, the EEG effect on the displacement of the r 
coordinate in the if) wave function, qAt/m, is on the order of 1.5 x 10 2 i m which 
within 100 ms is on the order of 1.5 x -3 m. If we assume the extent of Ar can be 
on the order of a synapse (40), then this spatial extent is on the order of about /im 
= 10 4 A(A= Angstrom = 10xl0~ 10 m). (Typical synaptic gaps are on the order of a 
few nm.) If this is correct, then the displacement of r by the A term is much larger 
than Ar. If the uncertainty principle is close to saturation, we can take Ap > h(2Ar) 
= 1.054 x 10~ 34 /(2 x 10~ 6 ) = 5 x 10 - 29 kg-m/s. This would make Ap about the same 
as p. Given this spread for most ions in a wave, it is reasonable to further investigate 
this "beam" of ions with respect to their entanglement. 



2. 5. Quantum Coherence of Ca 2+ Waves 



26 , the 



While Ca 2+ are observed to remain in waves for durations up to 500 ms 
example above invokes extremely long quantum coherence times of 100 ms, and even 
considerations on how long quantum coherence times may be achieved may not support 
these long times. In any case, it should be noted that A exerts strong quantum influences 
on r via its relative influence on p. However, there are reasons to consider effects that 
may promote long quantum coherence times. 

It is now understood that standard arguments, that quantum coherence cannot be 
maintained at high temperatures 49 , simply may not apply to many complex biological 
systems where other interactions may take precedence 50 -54 . Quantum coherence in 



potassium ion channels has been proposed 55 . 

Here, waves of free Ca 2+ ions (26], the ions being synchronized into coherent waves 

may introduce 



56 



by phase coordination in this columnar coherent firing context 
pulsed-dynamical decoupling, a generalization of the quantum Zeno effect (QZE) and 
"bang-bang" (BB) decoupling, of ions from their environment, promoting long coherence 
times 



57-59], as the system receives n "kicks" during time t, 
U n (p,t) = [U k U(p,t/n)} n 



(12) 



where the kicks may include interactions with other quantum systems, e.g., 
other Ca 2+ ions in the same wave. Such mechanisms for maintaining coherence are 
currently investigated in the context of quantum computation 60 61 . Distinguishable 



particles, sometimes even if previously uncorrelated, can exhibit quantum coherence 
and entanglement via collisions 62,63 , e.g., such as collisions via coulomb interactions 
among ions being synchronously influenced by A in Ca 2+ waves. 
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In the context of of quantum Ca waves interacting with A, while it is 
straightforward to model many-body Hamiltonians/Lagrangians and wave functions, 
this issue can only be resolved by experimental verification, e.g., to ascertain the degree 
of quantum coherence among ions in a Ca 2+ wave. There is as yet no experimental 
evidence as to how this long-time coherence is developed by Ca 2+ waves in vivo. 



3. Coupled SMNI-VP Ca 2+ - Waves 
3.1. SMNI Dipoles 

A dipole model for collective minicolumnar oscillatory currents is considered, 
corresponding to top-down signaling, flowing in ensembles of axons, not for individual 
neurons. The top-down signal is claimed to cause relevant effects on the surrounding 
milieu, but is not appropriate outside these surfaces due to strong attenuation of 
electrical activity. However, the vector potentials produced by these dipoles due to 
axonal discharges do survive far from the axons, and this can lead to important effects 



at the molecular scale, e.g., in the environment of ions 36 64 . 

The SMNI columnar probability distributions, derived from statistical aggregation 
of synaptic and neuronal interactions among minicolumns and macrocolumns, have 
established credibility at columnar scales by detailed calculations of properties of 
STM. Under conditions enhancing multiple attractors, detailed in SMNI with a 
"centering mechanism" (CM) effected by changes in background synaptic activity, 
multiple columnar collective firing states are developed. These minicolumns are the 
entities which the above dipole moment is modeling. The Lagrangian of the SMNI 
distributions, although possessing multivariate nonlinear means and covariance, have 
functional forms similar to arguments of firing distributions of individual neurons, so 
that the description of the columnar dipole above is a model faithful to the standard 
derivation of a vector potential from an oscillating electric dipole. 

Note that this is not necessarily the only or most popular description of 
electromagnetic influences in neocortex, which often describes dendritic presynaptic 



activity as inducing large scale EEG 65 , or axonal firings directly affecting astrocyte 
processes 66 . This work is only and specifically concerned with electromagnetic 
fields in collective axonal firings, directly associated with columnar STM phenomena in 
SMNI calculations, which create vector potentials influencing ion momenta just outside 
minicolumnar structures. 



3.2. SMNI Lagrangian 

A very short summary of the relevant SMNI Lagrangian in terms of its scaled synaptic 
parameters enables an explicit presentation of coupling the SMNI-VP with Ca 2+ waves. 

Care was taken in the first derivations of SMNI to properly process time-dependent 
and nonlinear multivariate drifts and diffusions. E.g., in the mid-point (Stratonovich or 
Feynman) representation, all Reimannian contributions were calculated and numerically 
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estimated for neocortex, as the nononlinear multivariate diffusions present a curved 



space 11,12 . A derivation of the underlying mathematical physics has been in some 



specialized text books for some time 67 , and a compact derivation has been given in 



several papers 68 70 



The SMNI Lagrangian, L, in the prepoint (Ito) representation was derived as 
L = ^(2iV)- 1 (M G - g G )g GG ,{M G ' - g G ')/{2Nr) - V 

G,G' 

g G = -t-\M g + N G tanhF G ) 
9 GG ' = (9GG')- 1 = S G 'r-'N G se C h 2 F G 

g = det(g GG ,) (13) 

where G = {E, 1} represents excitatory E and inhibitory I processes, the aggregated 
relaxation time r is on the order of 10 ms, N = N E + N 1 , and N E = 80, N 1 = 30 has 
been used the number of E and I neurons in a minicolumn, with twice these numbers for 
visual cortex, in previous calculations. V are derived mesocolumnar nearest-neighbor 
(NN) interactions among minicolumns within macrocolumns. 
The threshold factor F G is derived as 

v G + iA E ' 



F G = V 

^(( 7r /2)[(,;g,) 2 + (0gO 2 ](^ + ^')) V2 

V G = V G - a%,v%N G ' - ^A%,v%M G ' 

U XE> = - a f lV E N tE' _ l _ A \E v E iM XE> 

5 G = a G N G ' + \a g M g ' 

8& = a XE N XE> + l _ A tE M %E> 



r -E> 



a G = \a g , + B G , a% E = l -A^, + B% E (14) 

where Aq, and Bq, are minicolumnar-averaged inter-neuronal synaptic efficacies, v G , and 
(f>Qi are averaged means and variances of contributions to neuronal electric polarizations. 
M G ' and N G ' in F are afferent macrocolumnar firings, scaled to efferent minicolumnar 
firings by N/N* w 10 -3 , where A^* is the number of neurons in a macrocolumn, 
about 10 5 . Similarly, Aq' and Bq have been scaled by N * /N ~ 10 3 to keep F G 
invariant. Other values taken are consistent with experimental data, e.g., V G = 10 mV, 
v% = 0.1 mV, <$>% = 0.03 1 / 2 mV. 

The numerator of F G contains post-synaptic parameters, and that the denominator 
of F G contains pre-synaptic parameters, a result that drops out of the derivation of the 
mesoscopic derivation from the statistics of synaptic and neuronal interactions in and 
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across minicolumns. Afferent contributions from N^ E long-ranged excitatory fibers, 
e.g., cortico-cortical neurons, are included, where N* E might be on the order of 10% 
of N*: Of the approximately 10 10 to 10 11 neocortical neurons, estimates of the number 
of pyramidal cells range from 2/3 up to 4/5 [71]. Nearly every pyramidal cell has an 
axon branch that makes a cortico-cortical connection; i.e., the number of cortico-cortical 
fibers is of the order 10 10 . This development is used in the SMNI description of scalp 
EEG across regions. 



3.2.1. Euler-Lagrange Equations The Lagrangian components and Euler-Lagrange 
(EL) equations are essentially the counterpart to classical dynamics, 

3 2 L 



Mass = g G G' 



Momentum = IT G 



d(dM G /dt)d{dM G '/dt) 
dL 



Force 



d(dM G /dt) 
dL 



dM G 



dL d dL 

F-ma = 0: 5L = = ^ - ^^c/q^ ( 15 ) 

Concepts like momentum, force, inertia, etc., are so ingrained into our culture, that 
we apply them to many stochastic systems, like weather, financial markets, etc., often 
without giving much thought to how these concepts might be precisely identified. For a 
large class of stochastic systems, even including nonlinear nonequilibrium multivariate 
Gaussian-Markovian systems, like SMNI, the above formulation is precise. The 
Momentum defined above are also used as Canonical Momentum Indicators (CMI) in 
several studies that demonstrated its superiority over simple statistical correlations as 
they take into account some physical properties of the systems studied. 



3.3. Coupling Ca 2+ -Waves with SMNI Lagrangian 

The SMNI approach is a bottom-up mesoscopic aggregation from microscopic synaptic 
to columnar scales, and then scaled to relatively macroscopic regional scales of neocortex, 
which has been further merged with larger non-invasive EEG scales — all at scales 
much coarser than molecular scales. Here it is calculated how an SMNI vector potential 
(SMNI- VP) constructed from magnetic fields induced by neuronal electrical firings, at 
thresholds of collective minicolumnar activity with laminar specification, can give rise to 
causal top-down mechanisms that effect molecular excitatory and inhibitory processes 
in STM and LTM. 

While many studies have examined the influences of changes in Ca 2+ distributions 
on large-scale EEG 72 , future work will examine the complimentary effects on Ca 2+ 



ions at a given neuron site from EEG-induced magnetic fields arising from other neuron 
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sites. Here, sufficient calculations claim the importance of macroscopic EEG A, arising 
from microscopic synchronous neural activity, on molecular momenta p in Ca 2+ ions. 
The time dependence of Ca 2+ wave momenta may be calculated with rate-equations 



73 as a Hodgkin-Huxley model 74 , including contributions from astrocytes in the 



vicinity of synapses 75 , 76 . In this study, the resulting flow of Ca 2+ wave momenta will 
be further determined by its interactions in II, the canonical momenta which includes 
A. 



One influence of Ca 2+ likely is regulating synaptic interactions 77 . The SMNI 
Lagrangian explicitly describes where the Ca 2+ affect the columnar-averaged synaptic 
parameters Aq,, B g ,, A^, B^j . In this context Ca 2+ wave activity can affect the A and 
B synaptic parameters in these equations, while the A EEG fields affect the Ca 2+ waves. 

In SMNI papers, the CM is invoked by fine-tuning B parameters to bring maximum 
multiple minima in firing space M, by adjusting background Bq, to set v G = 
when M G = 0, similar to the control of spontaneous synaptic background observed 



during selective attention 78 . The B parameters are a logical first choice to include 
influences from columnar Ca 2+ activities. These minima tend to lie along a line in 
a trough in M space, essentially A E M E — AfM 1 « 0, noting that in F 1 I — I 
connectivity is experimentally observed to be very small relative to other pairings, so 
that (A I E M E — AjM ) is typically small only for small M E . This model gives rise to all 
the successful SMNI calculations describing various STM phenomena. 

This trough also supported previous SMNI work fitted to EEG data 16 , developing 
a scaled macrocolumnar electric potential at scalp region v derived with first and 
second moments of the SMNI Lagrangian, the argument of the associated distribution 
P u , at each (interconnected) region u, 

P v [$ v (t)\Mt - At)] = {27T J At)1/2 exp(-L„At) 
Lu= 2^ ( ^" m)2 

m =< - <p >= a < M E > +b < M 1 >= ag E + bg 1 

a 2 =< ($„ - <pf > - < $„ - > 2 = a 2 g EE + b 2 g n (16) 

in terms of M G -space drifts g G , diffusions g GG \ and an averaged reference (p. 

The same process supports the similar parameterization of A in these studies, i.e., 

A = cM E r + alM E i (17) 

where c and d are scaled to the order of 10 4 pA, as discussed above. This results in 
a Lagrangian L for the combined EEG-Ca 2+ system, e.g., considering |A| as primarily 
perpendicular to the scalp, 

1 



r(|A|„-m' 



'\2 



2a' 2 

m! =< \A\ U — <fi' >= c < M E > +d < M 1 >= cg E + dg 1 
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a' 2 =< (\A\ V - </>') 2 >-< \A\ U - (f)' > 2 = c 2 g EE + d 2 g n (18) 

where now m! and a' are nonlinear functions of A via the B synaptic parameters 
discussed further below. This Lagrangian is the argument of the exponential defining the 
conditional probability density for developing from a state at time t — 1 to time t. The 
variational principle obeyed by this Lagrangian permits optimization of parameters to 
find most likely states that best fit EEG data, i.e., including macrocolumnar parameters 



within regions, long-ranged connectivity and time delays across regions 16 



3.4- Experimental Verification 

The duration of a Ca 2+ wave can be on the order of 500 ms, so that the momenta of 
such ions can be importantly influenced during relatively long EEG events like N100 
and P300 potentials, reflecting latencies on the order of 100 ms and 300 ms, common 



in selective attention tasks which span these events 45 . Similar to procedures used 



in previous SMNI fits to EEG data 16 , 79 , here the influence of Ca 2+ waves may be 



tested by parameterizing the B synaptic parameters to include their influence in data 
sets where subjects have had simultaneous recording of scalp EEG and samplings of 
Ca 2+ wave activity at synaptic scales. These parameters are then fit to a portion of the 
EEG data, the in-sample set. The trained parameters can be be used in out of sample 
EEG data, to test if the included Ca 2+ wave activity correlates with the observed Ca 2+ 
wave data. 

The interaction of A and Ca 2+ waves can be detailed using SMNI-scaled synaptic 
parameters which include a term dependent on A, with coefficients measuring the 
convergence of synaptic interactions from many local minicolumnar and regional long- 
ranged fibers. The waves depend on aggregates of their II = p + gA dynamics. E.g., 
this can be modeled Taylor expansion in |A|, 

B% ->• B% + \A\B% , B|f = B %E + \A\B' tE (19) 

Eventually, the functional form of these dynamics should be established by models fit 
to molecular dynamics data, but for now at least their parameterized influences can 
be included. Since $ is experimentally measured, not A, but both are developed by 
currents I, when fitting to EEG data, it is reasonable to consider A as proportional to $ 
with a simple scaling factor, and now the additional parameterization of Bq, and B^ are 
to be included to modify previous work. To handle the otherwise recursive calculation 
of |A| multiplying B'q, and B'^j , here |A| is saved as a multiple of \g GG \r from just- 
previous data points, to be used in current time in the cost function calculation. A 
reasonable constraint is imposed that the inclusion of the B' terms not exceed the value 
of the B terms, e.g., limiting the influence of B' to at most doubling the background 
noise. The data used for this study, is spaced about 3.6 ms (< r) between 150-400 ms 



after presentation of stimuli 16 79 



The values of averaged synaptic parameters used in the 1980's SMNI papers were 
taken from experimental papers. Without any fitting of these parameters to other 
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data, SMNI detailed STM phenomena, e.g., as mentioned previously, capacity (auditory 
7 ± 2 and visual 4 ± 2), duration, stability , primacy versus recency rule, Hick's law - 
and other properties of neocortex by scaling up to macrocolumns across regions to 
fit EEG data. The S's terms, previous to the present B' inclusion, were consistent 
with these STM observations. Therefore, in this study, the B' terms were constrained 
to add no more than their B counterparts. Furthermore, since the data being fit is 
within the duration of P300 EEG waves, the inclusion of the time-dependent B' terms, 
i.e., including time-dependent modeled |A|, required a "dynamic centering mechanism" 
(DCM) to model regular access to maximum memory states, consistent with the early 
SMNI studies. Another study mentioned below will simulate the contribution of Ca 2+ 
waves via tripartite synaptic interactions, to determine if the changes implemented are 
reasonable assumptions. 

Another experimental test at the classical molecular scale to verify the influence of 
A, can be made considering that if the current lies along z, then A only has components 
along z, and 

II = p x x + py-y + (p x + qA z )z (20) 

The influence of time-dependent Ca 2+ waves is introduced in the post-synaptic and 
pre-synaptic SMNI parameters, which here also are time-dependent as functions of 
changing Ca 2+ ions. Such parameters are present at neuronal scales and are included 
in microscopic ordinary differential equation calculations. However, as in the original 
development of SMNI, these parameters are developed to mescolumnar scales, and the 
prediction here is that there is a predominance of Ca 2+ waves in directions closely aligned 
to the direction perpendicular to neocortical laminae during strong collective EEG. 

4. Current Projects 

The author is Principal Investigator (PI) of the National Science Foundation (NSF) 
collaborative resource The Extreme Science and Engineering Discovery Environment 
(XSEDE.org), working on this project, "Electroencephalographic field influence on 
calcium momentum waves". Some similar prior work was accomplished in 1994-1995 
when the author was PI of an NSF/PSC Grant, "Porting Adaptive Simulated Annealing 
and Path Integral Calculations to the Cray; Parallelizing ASA and PATHINT Project" 
(PAPP). 

Several sub-projects are being developed, using codes that permit some 
parallelization. 

4.1. EEG Data Fits 

EEG data is fit to SMNI, using data that the author has made public at 
kdd.ics.uci.edu/databases/eeg/ as described in| http: / /ingber.com/smni97_eeg-data.html} 
which is regularly used by other researchers. This project examines the influence of A 
on the B synaptic parameters in the SMNI Lagrangian. 



EEG influence on calcium waves 



15 



4-1.1. Data EEG spontaneous and evoked potential (EP) data from a multi-electrode 
array under a variety of conditions was collected at several centers in the United States, 
sponsored by the National Institute on Alcohol Abuse and Alcoholism (NIAAA) project 



80 -82 . This data set was used in earlier SMNI studies 16 , 79 . These experiments, 
performed on carefully selected sets of subjects, suggest a genetic predisposition to 
alcoholism that is strongly correlated to EEG AEP responses to patterned targets. 



The ASA code is used for fitting SMNI to this data 30 . These fits permit an 
estimate of the influence of |A| on the B synaptic parameters. It suffices for this 
study that the data used is clean raw EEG data, and that these SMNI, CMI, and ASA 
techniques can and should be used and tested on other sources of EEG data as well. 

Each set of results is presented with 6 figures, labeled as [{alcoholic | control}, 
{stimulus 1 I match | no-match}, subject, {potential | momenta}], abbreviated to {a 
I c}_{1 I m I n}_subject where match or no-match was performed for stimulus 2 after 



3.2 sec of a presentation of stimulus 1 80 82 . Data includes 10 trials of 69 epochs 



each between 150 and 400 msec after presentation. For each subject run, after fitting 
28 parameters with ASA, epoch by epoch averages are developed of the raw data and 
of the multivariate SMNI CMI. It was noted that much poorer fits were achieved when 



the CM 13 , 83 , driving multiple attractors into the physical firing regions bounded 
by M G < ±N G , was turned off and the denominators in F G were set to constants, 
confirming the importance of using the full SMNI model. All stimuli were presented 
for 300 msec. Note that the subject number also includes the {alcoholic | control} tag, 
but this tag was added just to aid sorting of files (as there are contribution from co2 
and co3 subjects). Each figure contains graphs superimposed for 6 electrode sites (out 
of 64 in the data) which have been modeled by SMNI using a circuitry of frontal sites 
(F3 and F4) feeding same-side sites: temporal sites (T7 and T8) with delay times of 1 
unit of data resolution, parietal sites (P7 and P8) with delay times of 2 units of data 
resolution. Temporal sites also feed same-side parietal sites with delay times of 1 unit of 
data resolution. Additionally, there are cross-side interactions, between temporal sites 
and between parietal sites with delay times of 1 unit of data resolution. Odd-numbered 
(even-numbered) sites refer to the left (right) brain. 

4-1-2. Optimization This optimization used ASA, 2013 version 28.15, tuned to give 
reasonable performance by examining intermediate results of several sample runs in 



detail. See the ASA code for a discussion of ASA OPTIONS and Tuning 30; . 

For both A with DCM and no-A models with CM, ASA was used for 60 data sets in 
{a_n, a_m, a_n, c_l, c_m, c_n} of 10 subjects. Each of these data sets had 4-6 parameters 
for each SMNI electrode-site model in {F3, F4, T7, T8, P7, P8}, i.e., 34 parameters 
(28 parameters when Bl terms are dropped) for each of the optimization runs, to be fit 
to over 400 pieces of potential data. This again is the same procedure used in previous 



papers with this data 16 79 



The ranges of the parameters were decided as follows. The ranges of the strength of 
the long-range connectivities d u were from to 1. The ranges of the {a, b, 0} parameters 
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were decided by using minimum and maximum values of M G and M^ G firings to keep 
the potential variable within the minimum and maximum values of the experimentally 
measured potential at each electrode site. 

It was found that typically within several thousand generated states, the global 
minimum was approached within at least one or two significant figures of the effective 
Lagrangian (including the prefactor). This estimate was based on fits achieved with 
2,000,000 generated states per run, after which ASA used its supplementary simplex 
code for additional local fits to sometimes get tighter fits. Each run generates 10 
files, of type asa-out, asa-usr-out, train-run-out, test-run-out, train-data, test-data, 
train-test-graph.ps, train-test-graph.eps, train-test-graph.pdf, and train-test-graph. png. 
Each ASA optimization took about 6 CPU-hrs for each of 120 runs on the XSEDE 
Trestles cluster, a cumulative CPU-month+ for 1200 files. An additional set of runs 
used 4,000,000 generated states per run, which improve a few graphs; this set is given 
here. 

Each complete set of runs were performed on the XSEDE Trestles machine, using 
four nodes, each node spawning 30 processors. Trestles is a dedicated XSEDE cluster 
designed by Appro and SDSC consisting of 324 compute nodes, running under Linux 
CentOS. Each compute node contains four sockets, each with a 8-core 2.4 GHz AMD 
Magny-Cours processor, for a total of 32 cores per node and 10,368 total cores for the 
system. Each node has 64 GB of DDR3 RAM, with a theoretical memory bandwidth of 
171 GB/s. The compute nodes are connected via QDR InfiniBand interconnect, fat tree 
topology, with each link capable of 8 GB/s (bidirectional). Trestles has a theoretical 
peak performance of 100 TFlop/s. 



4-1.3. Testing Data When the parameters of a theory of a physical system possess clear 
relationships to observed physical entities, and the theory fits experimental phenomenon 
while the parameters stay within experimentally determined ranges of these entities, 
then generally it is conceded that the theory and its parameters have passed a reasonable 
test. It is argued that this is the case for SMNI and its parameters, and this approach 
sufficed for the first study of the present data ll6j, just as SMNI also has been tested in 
previous papers. 

When a model of a physical system has a relatively phenomenological nature then 
often such a model is best tested by first "training" its parameters on one set of data, 
then seeing to what degree the same parameters can be used to match the model to out- 
of-sample "testing" data. For example, this is performed for the statistical mechanics 



of financial markets (SMFM) project, applied to trading models 84 85 



In the present project, there exists barely enough data to additionally test SMNI in 
this training versus testing methodology. That is, when first examining the data, it was 
decided to to try to find sets of data from at least 10 control and 10 alcoholic subjects, 
each set containing 10 runs for each of the 3 experimental paradigms, as reported in a 



previous paper 16 . When reviewing this data, e.g., for the example of the one alcoholic 



and the one control subject which were illustrated in graphs in that previous paper, 
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it was determined that there exists 10 additional sets of data for each subject for each 
paradigm, except for the c_n case of the no-match paradigm for the control subject where 
only 5 additional out-of-sample runs exist. For this latter case, to keep the number of 
runs sampled consistent across all sets of data, e.g., to keep the relative amplitudes of 
fluctuations reasonably meaningful, 5 runs of the previous testing set were joined with 
the 5 runs of the present training set to fill out the data sets required for this study. 



If.. 1.4- Graphical Results This procedure was followed both when including A terms 
in the synaptic background parameters in the SMNI model, and when excluding these 
terms (the previously published SMNI model). 

With the A terms, for both Training and Testing, Figs. 1-3 compares the CMI 
to raw data for an alcoholic subject for the a_l, a_m and a_n paradigms, for both the 
training and testing data; Figs. 4-6 gives similar comparisons for a control subject for 
the c_l, c_m and c_n paradigms. Without the A terms, for both Training and Testing, 
Figs. 7-9 compares the CMI to raw data for an alcoholic subject for the a_l, a_m and aji 
paradigms, for both the training and testing data; Figs. 10-12 gives similar comparisons 
for a control subject for the c_l, c_m and c_n paradigms. The CMI have been truncated 
to lie with [-0.5, 0.5] to keep the scale exposing their structure; an occasional small 
volatility can cause a point to jump out of this scale. To facilitate comparisons, the 
representative subjects presented here were selected to be the same as those presented 



in original studies 16 79 



The SMNI CMI clearly give better signal to noise resolution than the raw data, 
especially comparing the significant matching tasks between the control and the alcoholic 
groups, e.g., the c_m and a_m paradigms, in both the training and testing cases. The 
CMI can be processed further as is the raw data, and also used to calculate "energy" 
and "information/entropy" densities. 

Similar results are seen for other 10 control and 10 alcoholic subjects. For each set 
of |A| and no-|A| runs, there are 360 files of output data, 240 files of 4 graphs each (in 
ps, eps, png, and pdf formats), 80 files of parameter tables, and 13 files of summary 
statistics. 

After the above training-testing methodology is applied to more subjects, with 
additional variants of A models, it will then be possible to perform additional statistical 
analyses to seek more abbreviated measures of differences between alcoholic and control 
groups across the 3 experimental paradigms. 
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Figure 1. FIG. 1. With the A model, the initial-stimulus a_l paradigm for alcoholic 
subject co2a0000364, plots are given of activities under 6 electrodes of the CMI in the 
lower figures, and of the electric potential in the upper figures. 
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Figure 2. FIG. 2. With the A model, the match second-stimulus a_m paradigm for 
alcoholic subject co2a0000364, plots are given of activities under 6 electrodes of the 
CMI in the lower figures, and of the electric potential in the upper figures. 
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Figure 3. FIG. 3. With the A model, the no-match second-stimulus a_n paradigm 
for alcoholic subject co2a0000364, plots are given of activities under 6 electrodes of the 
CMI in the lower figures, and of the electric potential in the upper figures. 
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Figure 4. FIG. 4. With the A model, the initial-stimulus c_l paradigm for control 
subject co2c0000337, plots are given of activities under 6 electrodes of the CMI in the 
lower figures, and of the electric potential in the upper figures. 
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Figure 5. FIG. 5. With the A model, the match second-stimulus c_m paradigm for 
control subject co2c0000337, plots are given of activities under 6 electrodes of the CMI 
in the lower figures, and of the electric potential in the upper figures. 
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Figure 6. FIG. 6. With the A model, the no-match second-stimulus c_n paradigm 
for control subject co2c0000337, plots are given of activities under 6 electrodes of the 
CMI in the lower figures, and of the electric potential in the upper figures. 
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Figure 7. FIG. 7. With the no-A model, the initial-stimulus a_l paradigm for 
alcoholic subject co2a0000364, plots are given of activities under 6 electrodes of the 
CMI in the lower figures, and of the electric potential in the upper figures. 
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Figure 8. FIG. 8. With the no-A model, the match second-stimulus a_m paradigm 
for alcoholic subject co2a0000364, plots are given of activities under 6 electrodes of the 
CMI in the lower figures, and of the electric potential in the upper figures. 
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Figure 9. FIG. 9. With the no-A model, the no-match second-stimulus a_n paradigm 
for alcoholic subject co2a0000364, plots are given of activities under 6 electrodes of the 
CMI in the lower figures, and of the electric potential in the upper figures. 
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Figure 10. FIG. 10. With the no-A model, the initial-stimulus c_l paradigm for 
control subject co2c0000337, plots are given of activities under 6 electrodes of the CMI 
in the lower figures, and of the electric potential in the upper figures. 
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Figure 11. FIG. 11. With the no- A model, the match second-stimulus c_m paradigm 
for control subject co2c0000337, plots are given of activities under 6 electrodes of the 
CMI in the lower figures, and of the electric potential in the upper figures. 
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Figure 12. FIG. 12. With the no-A model, the no-match second-stimulus 
c_n paradigm for control subject co2c0000337, plots are given of activities under 6 
electrodes of the CMI in the lower figures, and of the electric potential in the upper 
figures. 
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4-1.5. Discussion While the CMI show marked differences between the A and no A 
runs, the degrees of fitting, measured by the values of the cost functions L, are close, 
but noA does seem marginally better (with 6 fewer parameters as well). No definite 
conclusion can be drawn for these runs regarding the superiority of either model. This 
likely is due to the constraint imposed on the A model, that the inclusion of the B' 
terms not exceed the value of the B terms, e.g., limiting the influence of B' to at most 
doubling the background noise; this saturation happened quite often. Other A models 
are being tested. 



4-2. Contribution of Ca 2+ Waves 

or some model 



The NEURON code 86 or some model 87 is used to simulate contributions of Ca 2+ 
waves from "internal stores" as discussed above. This permits an estimate of the 
importance of tripartite Ca 2+ waves to the the B synaptic parameters as derived in 
SMNI papers 11 12 . In those 1982 and 1983 papers, derivations based on different 



distributions for the number of quanta q released across synaptic gaps, gave the same 
means for the distribution of q, represented by a% and a^j. The a's are a sum of l/2A's 
due to neuronal firings plus background B's due to background noise. It is reasonable 
to consider that Ca 2+ waves from tripartite interactions contribute to the .B's. 



4-3. Quantum Coherence Simulations 



We are researching the calculation of extended quantum coherence of Ca using 
PATHINT 32 or PATHTREE 34 code used in previous publications by the author. 
The use of these codes for path-integral calculations, in contrast to Monte Carlo 
codes, permits a time step-wise propagation of quite general time-dependent nonlinear 
multivariate propagators, during which new events may enter the propagation, e.g., 
simulating BB decoupling from interacting ions in Ca 2+ waves, promoting long coherence 
times. This is quite a familiar situation in financial derivatives like options, and these 
codes were used successfully in that discipline, e.g., much better than using Monte Carlo 



calculations of path integrals 33 , 34 . However, here the propagator may not be easily 
defined (we have not found any good example yet), and it likely will in complex (x + iy) 
space which makes numerical details quite harder. 

Quite apart from obvious possible extrapolations to the realm of interactions 
between macroscopic processes involved in consciousness and the quantum scales 
of neocortical processes, the issues in this project are within currently accessible 
experimental and theoretical physics per se. 



5. Conclusion 



A model has been developed to calculate and experimentally test the coupling of 
molecular scales of Ca 2+ wave dynamics with A fields developed at macroscopic regional 
scales measured by coherent neuronal firing activity measured by scalp EEG. 
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For several decades biological and biophysical research into neocortical information 
processing has explained neocortical interactions as specific bottom-up molecular and 



smaller-scale processes 88 . It is clear that most molecular approaches consider it 
inevitable that their approaches at molecular and possibly even quantum scales will 
yet prove to be causal explanations of relatively macroscopic phenomena. 

This study crosses molecular, microscopic (synaptic and neuronal), mesoscopic 
(minicolumns and macrocolumns), and macroscopic regional scales. Over the past 
three decades, with regard to STM and LTM phenomena, which themselves are likely 
components of other phenomena like attention and consciousness, the SMNI approach 



has yielded specific details of STM not present in molecular approaches 17 . The 



SMNI calculations detail information processing capable of neocortex using patterns of 
columnar firings, e.g., as observed in scalp EEG 19 , which give rise to a SMNI vector 
potential A that influences the molecular Ca 2+ momentum p, and thereby synaptic 
interactions. Explicit Lagrangians have been given, serving as cost /objective functions 
that can be fit to EEG data, as similarly performed in previous SMNI papers. 

This study predicts a predominance of Ca 2+ waves in directions closely aligned to 
the direction perpendicular to neocortical laminae (A is in the same direction as the 
current flow, typically across laminae, albeit they are convoluted), especially during 
strong collective EEG (e.g., strong enough to be measured on the scalp, such as during 
selective attention tasks). Since the spatial scales of Ca 2+ wave and macro-EEG are 
quite disparate, an experimenter would have to be able to correlate both scales in time 
scales on the order of tens of milliseconds. 

The basic premise of this study is robust against much theoretical modeling, 
as experimental data is used wherever possible for both Ca 2+ ions and for large- 
scale electromagnetic activity. The theoretical construct of the canonical momentum 
II = p + qA is firmly entrenched in classical and quantum mechanics. Calculations 
demonstrate that macroscopic EEG A can be quite influential on the momentum p of 
Ca 2+ ions, at scales of both classical and quantum physics. 

A single Ca 2+ ion can have a momentum appreciably altered in the presence of 
macrocolumnar EEG firings, and this effect is magnified when many ions in a wave 
are similarly affected. Therefore, large-scale top-down neocortical processing giving rise 
to measurable scalp EEG can directly influence molecular-scale bottom-up processes. 
This suggests that, instead of the common assumption that Ca 2+ waves contribute to 
neuronal activity, they may in fact at times be caused by the influence of A of larger- 
scale EEG. Such a top-down effect awaits forensic in vivo experimental verification, 
requiring appreciating the necessity and due diligence of including true multiple-scale 
interactions across orders of magnitude in the complex neocortical environment. 
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